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Abstract Understanding the influence of macromolecular crowding and nanoparticles on tlie 
formation of in-register /3-sheets, the primary structural component of amyloid fibrils, is a first 
step towards describing in vivo protein aggregation and interactions between synthetic materials 
and proteins. Using all atom molecular simulations in implicit solvent we illustrate the effects of 
nanoparticle size, shape, and volume fraction on oligomer formation of an amyloidogenic peptide 
from the transthyretin protein. Surprisingly, we find that inert spherical crowding particles 
destabilize in-register /3-sheets formed by dimers while stabilizing /3-sheets comprised of trimers 
and tetramers. As the radius of the nanoparticle increases crowding effects decrease, implying 
smaller crowding particles have the largest influence on the earliest amyloid species. We explain 
these results using a theory based on the depletion effect. Finally, we show that spherocylindrical 
crowders destabilize the ordered /3-sheet dimer to a greater extent than spherical crowders, which 
underscores the influence of nanoparticle shape on protein aggregation. 
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Molecular crowding can have a profound effect on virtually all biological processes such as 
protein folding, viral capsid assembly and protein aggregation jl-10|. It is estimated that 20 



to 30% of a typical cell's volume is occupied by DNA, protein, lipids and sugars [3|, [ill, Il2 |. 
Such a crowded environment restricts the conformations explored by biological macromolecules, 
thus affecting the balance of thermodynamic forces that help regulate cellular processes. More 
recently, it has also been realized that synthetic nanoparticles (NPs) such as quantum dots, 
carbon nanotubes, and gold and other colloidal particles also affect the stability and function of 



proteins 
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15| . Understanding nanoparticle-protein interactions are important in our ability 



16|. 



to use NPs for drug delivery and controlling environmental toxicity 

In the context of aggregation of amyloidogenic peptides, which is the focus of the present 
study, NPs can greatly influence the stabilities of the molecular species that accumulate along the 
routes to fibril formation. Amyloid fibrils are experimentally characterized as cross-/? structures 



that are rich in /3-sheet content 
diseases 
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17H19|| . Formation of amyloid is linked to a number of human 



20| . and it appears that almost any polypeptide chain can form amyloid under 



appropriate solution conditions [isl. Even small peptides |2l|, such as the fragment comprised 
of residues 105-115 from the Transthyretin (TTR) protein, form amyloid fibrils {22!, 23|. TTR 
ibrils are associated with senile systemic amyloidosis and familial amyloid polyneuropathy I 
2Q|. 



I24j-l27l| suggest that monodis- 



General theoretical arguments based on the depletion effect 
perse spherical NPs should promote protein aggregation. The depletion effect, which results in 
non-specific entropy-induced attraction between proteins in the presence of non-adsorbing NPs, 



arises from the volume excluded to the NPs by the monomer subunits of aggregating proteins. 
Consider for example two proteins that are spherical with a radius . When two such pro- 
teins approach each other up to a distance on the order of the radius, Rc, of the crowding 
particles, the particles can no longer fit and are therefore expelled from the gap formed by the 
proteins. Thus, the crowding particles exert an osmotic pressure on the associating monomers 
resulting in a net inter-protein attraction 24j-|26| . As a result, inert crowdin g pa rticles promote 
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3Q|. 



protein-protein association, which is in accord with several experiments 

Although this qualitative argument is compelling, the molecular consequences of crowding 
effects on oligomer formation is unknown. In particular it is unclear whether crowding parti- 
cles stabilize ordered /3-sheet oligomers, which might facilitate amyloid formation, or stabilize 
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amorphous aggregates with httle /3-sheet content, which would hkely impede or delay amyloid 
formation. To address these unresolved issues we have carried out molecular simulations of 
systems composed of either two, three, or four TTR peptides to assess the effect of size of spher- 
ical NPs on the structural and energetic properties of TTR oligomers. We find that spherical 
NPs destabilize ordered dimer aggregates but stabilize trimer and tetramer ordered aggregates. 
Increasing the radius of the NP, at a fixed crowder volume fraction reduces the effect of 

molecular crowding due to an increase in the interstitial space between NPs. In contrast, in- 
creasing the aspect ratio of the NPs with a sphero cylinder shape leads to greater destabilization 
of the ordered dimer than a spherical crowding particle. 

In order to assess the effect of crowders on oligomer formation we first carried out separate 
constant temperature simulations on the monomer, dimer (denoted {TTR}2), trimer {{TTR}^), 
and tetramer {{TTR}^) systems in the EEFl.l implicit solvent model 3l| at = 0, a temper- 
ature of 395 K, and at a protein concentration around 30 mM (SI Table 1). A high temperature 
was used to increase conformational sampling and yield converged results; we found that simu- 
lations carried out at temperatures less than 330 K were not converged. It should be stressed 
;hat for a complete understanding of crowding effects peptide concentration must also be varied 
4 1 as has been recently demonstrated in a complementary study using model systems I32I. For 
the systems containing two or more peptides we find that due to the high protein concentration 
the peptides are associated. Therefore, we cannot examine the effect of NPs on the association 
process of these peptides. However, as mentioned previously, the more crucial aspect from the 
perspective of amyloid formation is determining the effect of NPs on the different aggregated 
species (ordered versus disordered), an issue which these simulations can directly address. 

The TTR monomer is compact and devoid of persistent secondary structure at = 0. Its 
average radius of gyration is 7 A, whereas a fully extended TTR structure would have a value 
of around 12 A. Analysis of the secondary structure content shows that coil and turn structures 
are the most prevalent while helical and /3-strand structures are negligibly populated (< 2%). 
Even in the presence of crowding agents the strand content does not change appreciably. Thus, 
in isolation the TTR peptide is unstructured. 

The dimer {TTR} 2 has an appreciable population of both ordered and disordered aggre- 
gates. In the ordered state the peptides prefer parallel in-register /3-sheets over anti-parallel 
structures. The fraction of /3-strand content in {TTR}2 is 0.42 (Fig. [1^). The value of P/r, 
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FIG. 1: Effect of spherical crowding on {TTR}n as a function ^c- (A) and (B) are for the {TTR}2- (A) 
Secondary structure content (/3-strand (black), random coil (red), or turn(green)) versus ^c- Example 
structures from the ordered aggregate and disordered aggregate are displayed on the upper left and 
right respectively. (B) The fraction of in-register contacts (Pir, see Eq. 4 in the SI) versus ^c- (C) 
and (D) correspond to data from the {TTR}^ and {TTi?}4 systems. (C) The ^c-dependent changes 
in the /3-strand content for the trimer (solid black line with circles) and the tetramer (dashed red line 
with squares). (D) Same as (C) except the ordinate is the probability of finding an in-register aggregate 
(including both parallel and anti-parallel configurations). 



the probability that the two strands are in-register, when both are in a parallel arrangement 
is nearly 3-fold greater than the probability of forming in-register anti-parallel structures (Fig. 
Wp)- Expanded peptide conformations have several more in-register contacts on average, make 
more inter-peptide backbone hydrogen bonds, and are more likely to have /3-strand content 
than compact peptide conformations. For example, the radius of gyration (Rg) of the individual 
peptides in the ordered {TTR}2 ('ordered' being defined as having greater than eight in-register 
contacts out of a possible eleven) is 10.8 A while {Rg) of the peptide in the disordered {TTR}2 
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conformations (which have less than 5 in- register contacts) is 6.7 A. 

Independent simulations of the {TTR}^ and {TTR}^ at = reveal they also exist 
as either disordered aggregates or ordered /3-sheets. However, the probabilities of /3-sheet 
structure, which includes out-of-register /3-strands, for {TTR}^ and {TT_R}4 are 0.52 and 
0.44, respectively (Fig. [Th). The probability that {TTR}^ and {TTi?}4 form ordered in- 
register parallel and anti-parallel strand arrangements are 0.30 and 0.06, respectively (Fig. [TJl). 
The decrease in ordered in-register /3-sheets, 
relative to the {TTR}2 system, is due to an 
increase in the number of low energy out- 
of-register disordered conformations that are 
accessible as the number of peptides in the 
aggregate increases. Thus, the landscape of 
TTR oligomers has a number of distinct basins 
of attraction with conformationally heteroge- 
neous structures (see also [33]). 

We next examined the influence of inert 
spherical NPs with Rc = 3.5 A ($c > 0, 
Fig. 12^). For the {TTR}2 system the pres- 
ence of spherical NPs decreases the stability 
of the expanded (conformations with Rg > 10 
A) dimer structures monotonically as in- 
creases, while the stability of the collapsed 
structures (i.e., Rg < 7.5 A) increases (Fig. 

Wp)- As a consequence, the average solvent ac- 

, „ 1 , , , , , figuration of the NPs durine; the simulations. (B) 

cessibie surface area and the molecular volume 

r v 1 ^ • / 1 , The free energy profile as a function of the radius- 

of the dimer decreases as mcreases (data 

not shown) . In addition, the /3-strand content of-gyration of the monomers in the dimer at various 
(Fig. [1^) and the probability of finding poly- ^'^ values. 

merization competent in-register dimers {Pm) decreases with increasing (Fig. Wp)- Thus, 
relative to $c = 0, crowding destabilizes in-register /3-sheet aggregates causing the equilibrium 
to shift towards collapsed associated dimer structures at all values of (see the structures in 




FIG. 2: Effect of spherical crowding on the size of 
the {TTR}2 dimer. (A) Snapshot of a single con- 



6 



Fig. m- 

Surprisingly, we find the opposite result for the {TTR}^ and {TTi?}4, whose ordered struc- 
tures are stabilized as increases. The /3-strand content in {TTR}^ remains the same as 
increases, whereas there is a slight increase for {TTR}4 (Fig. In addition, Pj^ increases 
for the trimer and tetramer with increasing (Fig. [Hi). The observed variations in crowding- 
induced changes in the stabilities of the ordered /3-sheet aggregates may be relevant to in vivo 



amy 
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oia 
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formation as /3-sheets are the primary structural component found in mature fibrils 
. In particular, NPs may inhibit the ordered dimer structures but stabilize higher 



order oligomers. 

This non-monotonicity in the n-dependence of the stability of the ordered aggregate of 
{TTR}n {n = 2,3,4) at various can be explained using a theory based on depletion forces. 



Scaled Particle Theory 37H39j suggests that a factor infiuencing crowder effects is the free en- 



ergy cost of creating a void within the configuration of NPs to accommodate {TTi?}„ when 
$(7 7^ 0. The free energy is a function of volume of {TTR}nS volume, which increases as n 
becomes larger. Clearly, the probability of finding a volume to accommodate {TTi?}„ decreases 
as n increases. Depletion-induced changes in the free energy of the ordered aggregate relative 
to the disordered aggregate, denoted AGn, is therefore proportional to IniVoA — Vda), where 
VoA {Vda) is the volume excluded to the NP's center-of-mass (also known as the 'depletion vol- 
ume') by the ordered (disordered) aggregate. Peptide conformational entropy changes also likely 
contribute to the behavior of AG„ with n. However, it is difficult to accurately estimate these 
entropy changes in the ordered and disordered aggregates due to crowding. In addition to con- 
formational entropy differences, stability is also determined by favorable enthalpic inter-peptide 
interactions. Enthalpic stabilization occurs for both the ordered and disordered aggregates as n 
increases because of an increase in the number of inter-peptide contacts. The specific value of 
n which stabilizes the ordered form should depend on the size and sequence of the peptide. For 
{TTR}n it is only when n becomes larger than 3 that the ordered aggregate structures become 
more stable in the presence of spherical NPs (Figs. and[T]i). 

If depletion forces are a major factor in the non-monotonic behavior of AAG^ with n then 
we can account for the simulation results usiiig the changes in the volume excluded to the 
NPs. From the Asakura-Oosawa (AO) theory 24] and a microscopic formulation I25I], as well 



as several other approaches 



37H39|], it follows that the difference in free energy between species 



7 



i and j is AGi{n) = Piy^^^n) — Vj^{n)), where P = ksT^c/Vc, is the osmotic pressure 
under ideal solution conditions, ks is Boltzmann's constant, Vc is the molecular volume of a 
nanoparticle, and V^^{n) is the volume excluded to the spherical nanoparticle by species i. There 
are three relevant species in the {TTR}n system; the soluble non-associated monomers (SM), the 
disordered aggregate (DA), and the ordered aggregate (OA). We calculate V^^{n) by assuming 
that the soluble monomers and disordered aggregate are hard spheres whose radii depend on n. 
In this case 

V!f{n) = ^[(1.927iVr + Rc? - {l-927Nir] (1) 
V^An) = Y[(l-927(n ■ iV,J°-6 + Rcf - (1.927(n ■ KaD^ (2) 

where Naa is the number of amino acids in the peptide and the exponent of 0.6 is the Flory 
scaling exponent that characterizes the size of a protein in a good solvent 40|]. The shape of 
the ordered /3-sheet aggregate is taken to be a stacked rectangular parallelepiped, where each 
parallelepiped corresponds to one TTR peptide. For this species Vex is computed using 

Kx^(^) = {l + 2Rc){w + 2Rc){h + 2Rc)-lwh-{8Rl-^7rRl), (3) 

where I, w, and h correspond to the length, width and height of a /3-sheet made up of n peptides, 
respectively. Therefore, we define / = 3.2 ■ Naa A, corresponding to the length of an extend /3- 
strand made up of Naa amino-acids; w = 8 K which corresponds to the distance that the side 
chains stick out above and below the /3-sheet; h = A.S^n A, where 4.85 A corresponds to the 
experimentally determined distance between the neighboring strands in a /3-sheet. 

Using Eqs. [T]|3]it is straightforward to compute the Vex for each species at different $(7, and 
thereby estimate AGi{n) at a given ^c- We find that indeed for a range of Rc values this model 
displays non-monotonic behavior (Fig. [3^) consistent with the trends observed in the simula- 
tions; the ordered dimer is destabilized by spherical crowders, but as more peptides are added 
the NPs stabilize the ordered /3-sheet aggregates. The physical origin of this unusual behavior 
arises because disordered dimer aggregate excludes less volume than the ordered dimeric /3-sheet 
(Fig. [3b), but as additional peptides are added to the system the volume per peptide increases to 
a greater extent in the disordered aggregate than the ordered /3-sheet (the same observation was 



made in ref. 



4l|). As a consequence, the free energy of the system is minimized when the dimer 



is disordered whereas NP-induced ordered trimer and tetramer structures are more stable. Reca- 
pitulation of the qualitative trends observed in the simulation model suggests that the depletion 
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Disordered Aggregate 
(3-sheet 



effect explains the changes in nanoparticle- 
induced changes in stabihty as n varies. 

There number of other parameters be- 
sides $(7 that also influence oligomerization. 
What happens when the radius of the NPs 
increases? For a fixed peptide concentration 
([Cp]) and $(7, larger NPs exclude less vol- 
ume to the protein than smaller crowders be- 
cause the interstitial space between crowders 
increases with crowder size. It can be shown 
that the strength of the depletion force is pro- 
portional to Therefore, larger crowding 
particles should have a smaller effect in alter- 
ing the stability of {TTi?}„. Explicit simu- 
lations that we ran of {TTR]2 at $c = 0.15 FIG. 3: (A) Theoretical prediction of the stability 
using three spherical NP sizes {Rc = 3.5, 6 of the ordered /3-sheet aggregate relative to the dis- 
and 11 A) show that indeed as Rc increases ordered aggregate as a function of the number of 
larger crowders shift the equihbrium towards TTR peptides in the system. The solution condi- 
expanded structures just as observed in bulk tions in this theoretical treatment are T = 395 K, 
simulations (Fig. i]). However, in the simu- $^7=0.15, and i?c=30 A. (B) Volume excluded (i.e., 
lations we observe that the average /3-sheet the depletion volume) to the NPs center of mass of 
content does not exhibit an increase as the ^he soluble monomers (green triangles), disordered 
size of the NP increases and the probability ^gg^egate (red squares) and ordered /3-sheet aggre- 
(^ 0.06) of finding ordered in-register dimers g^^^ ^^^^^^ ^-^^^^^^ calculated from Eqs. m 
is also unchanged. Thus, increasing Rc of 

spherical NPs decreases the impact of crowding - a result consistent with other theoretical 



42|. 



and simulation modeling 

Another variable that controls the 
oligomer formation is the shape of 



extent or 
the crowding 



ordered 
particles. 
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Indeed, in vivo most naturally occurring 
molecular crowders are unlikely to be spheri- 
cal. Analysis of protein structures in the Pro- 



tein Data Bank 



43| show that even in the 



folded state proteins adopt anisotropic shapes. 
To explore the effect of NP shape on dimer 
stability we simulated the {TTR}2 system in 
the presence of spherocylindrical NPs in the 
isotropic phase whose aspect ratio i=L/ D + 1, 
where L is the spherocylinder length and D 
is its diameter) is 3.3. We find that for the 
same value (below $(7=0.1, which is be- 
low the isotropic to nematic transition point) 
spherocylindrical NPs cause a greater loss of 
/9-strand content in the dimer than spheri- 
cal NPs. Thus, spherocylindrical NPs desta- 




11 12 13 



R,(A) 
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FIG. 4: (A) Effect of spherical NP size on {TTR}2 
at = 0.15. Changes in the free energy profile as a 
bilize ordered oligomers. This suggests that function of Rg of the monomer (within the dimer) as 
anisotropic NPs may have a greater impact the radius Rc of the spherical NP is varied. (B) Av- 
on amyloid formation than spherical NPs, and erage /3-sheet content of the {TTR}2 system versus 
that crowder shape is another important vari- for spherical crowders (black line and circles), 
able in understanding amyloid formation in and spherocylindrical crowders with an aspect ra- 
vivo. To further explore realistic models of tio of 3.3 (blue hne and squares). The inset shows 
oligomerization in the cellular context it would a spherocyhnder used in the simulations, 
be important to consider a soup of proteins as 



was done recently to explore folding and diffusion in models of E. Coli. 44 1. 

Although the systems investigated in this study are caricatures of cellular crowding they 
illustrate the complexity of peptide aggregation under in vivo conditions. While generic argu- 
ments suggest that depletion forces should promote protein association, our work shows that 
the structures of the oligomers can be dramatically altered by the size, shape, and volume 
fraction of the nanoparticles. Our study shows that an interplay of a number of factors deter- 
mines the equilibrium between the ordered and disordered oligomeric structures. For spherical 



10 



NPs the variables that determine the depletion forces, and hence the strength and range of the 

Rg_ 

Rc 



entropically-induced inter peptide attraction, are q = ^ (where Rg is the size of the peptide 



monomer), [Cp] the peptide concentration, and ^c- 

This predicted complexity in the crowding-induced diagram of states of the amyloidogenic 
peptide is not surprising given that in the field of polymer physics it was discovered previously 
that the phase diagram of a mixture of hard sphere colloidal particles and self-avoiding polymers 
in an athermal solvent is complicated, and depends on polymer concentration, $(7 and q [26!]. 
For the TTR peptides and Rc values examined here q ranges in value from 0.6 up to 2 (Note that 
Rg ~7 A for the peptide monomer). It follows from our work that larger q values can promote 
extended peptide conformations with high /3-strand content in aggregates consisting of more 
than two peptides. Conversely, in the extreme case of g << 1 (in the so called colloid limit [26]) 
we predict that the influence of crowding effects will be decreased, and the equihbrium between 
ordered and disordered species may not be significantly perturbed relative to bulk ($^7=0). 

The results presented here show that understanding of NP effects on protein aggregation 
requires determination of a global phase diagram in terms of a number of variables. The 
non-monotonicity in NP-induced changes in stability should be contrasted with the effect of 
crowding on monomeric protein folding stability, which monotonically increases in the presence 
of non-adsorbing crowding particles [4]. Further insights into in vivo aggregation will require 
extension of the present work by taking crowder phase and polydispersity into account perhaps 



along the lines suggested recently 12j. In such a complex environment polydispersity and 



451 ] ■ which could 



concentration fluctuations could drive phase separation in the macromolecules 
add additional complications in the process of protein aggregation. 

Methods 



Models for peptide and solvent. We chose the peptide fragment (sequence Tyr-Thr-Ile- 
Ala-Ala-Leu-Leu-Ser-Pro-Tyr-Ser) from the Transthyretin protein (TTR) corresponding to 



residues 105-115. The structure of these 
(PDB ID IRVS) using sohd state NMR 



2 residues in an amyloid fibril has been determined 



22 



231 ]. We cap the peptide's N-terminus with an 



acetyl blocking group and its C-terminus with an amine blocking group {NH2). For the peptide 
we use an all-atom representation except for non-polar hydrogen atoms, which are omitted in 
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the calculations. The EEFl.l implicit solvent model 
in conjunction with the CHARMM19 force-field 31 |. 



46| is used to account for solvent effects 



The EEFl.l implicit solvent model. The parameterization of the EEFl.l 3l| force-field 
requires the use of specific cutoff distances for non-bonded interactions. The Lennard- Jones 
(LJ) and electrostatic interactions are truncated at 9 A, with a switch function applied to the 
LJ term starting at 7 A. The Lorentz-Berthelot mixing rules |43| are used to determine the 
undefined cx-values between the atomic centers in the peptide. A distance dependent dielectric 
constant is used for electrostatic interactions. 

The LJ interactions in the EFFl.l protein parameters are restricted to a distance less than 9 
A, which prevents us from studying crowders with Rq values greater than 4.5 A in the standard 
Charmm code. We modified the CHARMM code to allow the LJ cutoff to be dependent on the 
identity of the interacting atomic centers, thereby allowing us to use larger crowding particles. 
For protein-protein interactions nonbonded interactions are truncated at 9 A, but for crowder- 
crowder and crowder-protein interactions we use a cutoff of 2Rc + '2 A and -Rc + 2 A, respectively. 
These procedures allow us to maintain the EEFl.l nonbond requirements for protein-protein 
interactions while allowing for larger crowding particles to be simulated. 

Simulation details. Simulations were carried out in the NVT ensemble at 395 K using 
Langevin dynamics with a friction coefficient of 1 ps~^ . This high temperature allowed the 
the dimer simulations to reach equilibrium, whereas simulations at 330 K were found to not 
have converged on our simulated time scales (data not shown). The SHAKE algorithm was 
used to fix the bond lengths of covalently bonded hydrogen, allowing the use of a 2 fs time-step. 

Typically, we generated ten independent trajectories for each crowder size and (^c- At least 
half of the trajectories were started with peptide conformations taken from equilibrated bulk 
simulations at $(7 = 0.0. When possible the other five starting conformations were taken from 
solution conditions that were closest to those of interest. For example, when simulating crowders 
with Rq = 6 and 11 A, half of the initial protein conformations were from the equilibrium 
simulations of crowders with Rc = 3.5 A and $(7 = 0.15. 

A box length of 60.0 A to 80.0 A is used for all spherical NP conditions, which gives a peptide 
concentration in the range from 15 mM to 31 mM. At $0 = 0.05, 0.10, 0.15, 0.20 there are 60, 
120, 180 and 240 crowders in the periodic box, respectively. Because correlations in the pair 
distribution function between crowders do not persist for more than half the box length, finite 
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size effects are minimized (Fig. SI). A summary of tlie simulation details can be found in Table 
SI. 

Models for spherical particles. Spherical NPs are modeled as Lennard- Jones particles. The 
interaction between sites i and j on two distinct NPs is 

The choice of e = 10^^^ kcal/mol makes Eq. H] essentially repulsive, and the condition VLj{rij = 
2Rc) = kBT, allows us to solve for a. We consider three values for Rc = 3.5, 6, and 11 A, which 
leads to cr= 108.7, 186.3, and 341.6 A, respectively. 

Model for spherocylinders and crowder-protein interactions. We define all protein-crowder a- 
values such that when a crowding particle interaction site is {Rc + 1) A from a protein atomic 
center Vlj{Rc + 1) = ksT. The resulting a values are 139.7, 217.4 and 372.6 A for Rc = 3.5, 
6 and 11 A, respectively. 

The covalent bond between the tethered spheres that comprise the spherocylinder was mod- 
eled using 

Vnir) = KBl2{ro-rf (5) 

where Kb = 30 kcal/(mol-A^) and = 4 A is the equilibrium bond length. The rigidity of the 
spherocylinder is maintained using a bond angle potential 

Va{9) = ^{9o-ef (6) 

with 9o = 180° and Ka = 0.6092 kcal/ degree^. An important characteristic of a spherocylinder 
is the aspect ratio A. A value of A = 1 corresponds to a spherical crowder with diameter 
D = 2Rc. In our study all spherocylinders have Rc = 3.5 A, and Lc, which is proportional 
to the number of spheres (A^'^) in the spherocylinder. Lc = {Nb — 1)4.0 A. We used A = 3.3, 
corresponding to a spherocylinder made up of 5 beads. 
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